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Abstract 

A stochastic cellular automata (CA) model for pedestrian dy- 
namics is presented. Our goal is to simulate different types 
of pedestrian movement, from regular to panic. But here we 
emphasize regular situations which imply that pedestrians an- 
alyze environment and choose their route more carefully. And 
transition probabilities have to depict such effect. The poten- 
tials of floor fields and environment analysis are combined in 
the model obtained. People patience is included in the model. 
This makes simulation of pedestrians movement more real- 
istic. Some simulation results are presented and comparison 
with basic FF-model is made. 

INTRODUCTION 

Modelling of pedestrian dynamics is actual problem at 
present days. Different approaches from the social force 
model ([2 1 and references therein) based on differential equa- 
tions to stochastic CA models (f5] |4] |6l and references 
therein) are developed. They reproduce many collective prop- 
erties including lane formation, oscillations of the direction at 
bottlenecks, the so-called "faster-is-slower" effect. These are 
an important and remarkable basis for pedestrian modelling. 
But there are still things to be done in order to reproduce in- 
dividual pedestrian behavior more realistic and carefully. 

The model presented takes its inspiration from stochastic 
floor field (FF) CA model [5|. Here a static field is a map that 
pedestrian may use to orient in the space. Dynamic field is 
used to model herding behavior in panic situations. 

It's known that regular situations imply that pedestrians 
analyze environment and choose their route more carefully 
(see Q and reference therein). Pedestrians keep a certain dis- 
tance from other people and obstacles. The more hurried a 
pedestrian is and more tight crowd is this distance is smaller. 
We adopted a mathematical formalization of these points 
from H. 

Pedestrians minimize efforts to reach their destinations: 
feel strong aversion to taking detours or moving opposite to 
their desired direction. However, people normally choose the 
fastest rout but not the shortest. This means that opportunity 
to wait (to stay at present place) has to be realized in the 
model. (Models [4, 5 1 (and other CA models) imply that peo- 



ple can stay at present place if there is no space to move only.) 
We realize this point (people patience) in the algorithm. 

As well as it's necessary to take into account that some 
effects are more reside for certain regions. For instance, clog- 
ging situations are more pronounced in the nearest to an exit 
areas. This means that spatial adaptivity of correspondent 
model parameters to be introduced in the model. All these 
changes and additions extend basis FF model towards emo- 
tional aspect and improve and make flexible decision making 
process. By this reason model obtained was named as Intelli- 
gent FF model. 

INTELLIGENT FLOOR FIELD MODEL 
Space structure 

As usual for CA models the space (plane) is sampled into 
cells 40cm x 40cm (it's an average space occupied by a pedes- 
trian in a dense crowd [ 3 1) which can either be empty or oc- 
cupied by one pedestrian (particle) only. 

The von Neumann neighborhood is used. It implies that 
each particle can move to one of four its next-neighbor cells 
(i.j) or to stay at the present cell at each discrete time step 
t — ► t + 1, e.L, Vmax — !■ (Empirically the average velocity of 
a pedestrian is about 1.3m/s. So real time corresponding to 
one time step in the model is about 0.3s. ) Such movement 
is in accordance with certain transition probabilities that are 
explained below. 

Floor fields 

Static (S) and dynamic (D) floor fields are introduced and 
discussed in fl]|3]|2]. For each cell (z, j) values of Sjj and £),,■ 
are given. 

Static floor field S describes the shortest distance to an exit 
(or other destination point that depends on a task). It doesn't 
evolve with time and isn't changed by the presence of the 
particles. The value of Sij is set inversely proportional to the 
distance from the cell (z, j) to the exit. One can consider S as 
a map that pedestrian can use to move to the target point, e.g., 
exit. 

Dynamic floor field D is a virtual trace left by the pedes- 
trians similar to the pheromone in chemotaxis. It is used to 
model a " long-ranged" attractive interactions between the 
pedestrians, e.g., herding behavior that is observed in panic 
situations. Dynamic floor field is time dependent. In each 
time step each Dij decays with probability 8 and diffuses 



with probability a <G [0, 1] to one of its four neighboring 
cells. Decay and diffusion lead to broadening, dilution, and 
finally vanishing of the trace. At t = for all cells Djj = 0. 
Dij — > Dij + 1 just after particle left cell 

Field S and D works in such a way that transition proba- 
bility increases in the direction of higher fields S and D, i.e. 
motion in such direction is more feasible. 

Environment analysis 

To make next step pedestrian observes surroundings. Let 
r > be maximum distance (in cells) at which pedestrian can 
see around, by 

_ J 1, cell is occupied by a pedestrian; 
* lJ ' ~ y 0, cell (ij) is empty; 

denote the occupation number. 

If desired direction (cell) is occupied to minimize efforts to 
reach destination or to realize the fasters rout pedestrian has 
to a have an opportunity to stay at present cell or move to one 
of the rest unoccupied nearest cells. 

Keeping apart from other people and obstacles can be sim- 
ulated by decreasing probability for such direction (let a be 
a name of the direction). For present position of some pedes- 
trian let cell be a next-neighbor in the direction a. For 
cell we will calculate a term [4| 

A U = ( l --( E fml+r-rfj)), 

m=i (/=;') 

where r* — distance to a nearest obstacle (wall, column, etc., 
but not pedestrian) in the direction a starting from cell (i,j) 

( r tj < r ); L fmi — number of pedestrians that are in 

m=i (l=j) 

the direction a starting from cell up to the nearest ob- 
stacle (sum is over lines (m) or over columns (I) depending 
on direction a), all cells behind obstacle are considered as 
occupied. 

Update rules 

Update rules for CA are the following: 

1 . For each pedestrian the transition probability py to move 
to cell one of four next-neighbors is 

Pij = Norm- 1 [(l--( £ fmi + r-rtj)) 

m=i (/=;') 

exp(k s Sij)exp(k D Dij) exp(fcj)] , (1) 

where ks,kp — sensitivity parameters; exp(Aj) — inertia 
effect: k\ > for the direction of pedestrian's motion in 



the previous time step and kj = for other cells; normal- 
ization 

Norm = Y J Ai j e ksS "e kDD 'ie k ' , 

(iJ) 

where the sum is over all possible target cells. 

2. If Norm = then pedestrian stays at the present cell, oth- 
erwise pedestrian chooses randomly a target cell 
based on the transition probabilities determined by (1). 

3. If Norm ^ and (1 - ftj) = (i.e., cell (ij)* is oc- 
cupied) then pedestrian chooses randomly a target cell 
once again. Now target cell is chosen among the follow- 
ing candidates: rest next-neighbors available for moving 
(i.e., (1 — fij) 7^ 0) and the present cell. For the avail- 
able next-neighbor cell new transition probability 
is ^f'p and pedestrian can stay at present cell with 

dJ) ''' 

probability ^ ^ (where sum is over all available can- 

(>■,/) ''' 

didates at now). Obviously, if there is no available next- 
neighbors then particle stays at present cell. 

In contrast to [5] this step gives an opportunity for 
pedestrians not to move and wait when preferable direc- 
tion will free. 

4. Whenever two or more pedestrians have the same target 
cell, the movement of all involved pedestrians is denied 
with probability fy, i.e. all pedestrians remain at their 
old places [5]. One of the candidates moves to the de- 
sired cell with the probability 1 — fry. Pedestrian that 
is allowed to move has the largest probability among 
all candidates (in this case probability is a measure of 
pedestrians physical strength). The other probabilistic 
method [ 1 , 4, 5 1 can be used here as well. 

5. Pedestrians that are allowed to move perform their mo- 
tion to the target cell. D at the origin cell of each 
moving particle is increased by one D, ; — > D, 7 + 1 and 
therefor can take any non-negative integer value. 

These rules are applied to all particles at the same time, i.e., 
parallel update is used. 

Model parameters 

There are several parameters in the model obtained. Their 
values and physical meaning are presented below. 

• r > — maximal distance at which pedestrian can feel 
the surroundings. People avoid to walk close to obsta- 
cles and other people. Ay = if there is no free space 
to move, Ajj = 1 if direction is free, < Ay < 1 for any 
other intermediate situation. 



• k$ > — sensitivity parameter that can be interpreted 
as the knowledge of the shortest way to the destination 
point, or as a wish to move in a certain direction. k$ — 
means that pedestrian don't use information from the 
field S. The higher k$ is movement of the pedestrians is 
more directed. 

• ko > — sensitivity parameter that can be interpreted 
as a rate of herding behavior. It is known that people 
try to follow others particulary in panic situations (2). 
ko = means that pedestrian chooses a way of their own 
ignoring ways of others. The higher ko is the herding 
behavior of the pedestrians is more pronounced. 

• k\ > — parameter that determines the strength of iner- 
tia which suppresses quick changes of the direction. 

• 0<foj < 1: 

f^VjU. iffe^O, 

I maxS,,^' J ' ' 

m = \ « 

I /j, otherwise, 

where /j S [0, 1] — friction parameter that controls the 
resolution of conflicts in clogging situations, pjj works 
as some kind of local pressure between the pedestrians. 
The higher jjjj is pedestrians are more handicapped by 
others trying to reach the same target cell. Such situa- 
tions are natural and well pronounced for nearest to exit 
(destination point) space. For other areas it's not typical 

but it's possible. So to realize it and make simulation of 

s- 

individuals realistic the coefficient — is introduced 

max S,- ; 

i,j 

(in contrast with original FF model |5 1). 

• 8, a € [0, 1] — these constants control diffusion and de- 
cay of the dynamic floor field [5|. "It reflects the ran- 
domness of people's movement and the visible range of 
a person, respectively. If the room is full of smoke, then 
8 takes large value due to the reduced visibility. Through 
diffusion and decay the trace is broadened, diluted and 
vanishes after some time." 

DISCUSSION OF THE MODEL 

Model obtained is simple. For one time step there are 0(n) 
calculations if n pedestrians are involved. It gives advantage 
over continuous social-force model [2] where each time step 
0(n 2 ) interaction terms have to be evaluated. The discrete- 
ness of the model is advantage as well. It allows for a very ef- 
ficient implementation for large-scale computer simulations. 

It's shown [TJ O that original FF model reproduces va- 
riety of collective effects: clogging at large densities, lane 
formation in counterflow, oscillation in counterflow at bot- 
tlenecks, patterns at intersection, trail formation, "faster-is- 



slower" and "freezing-by-hearting" (in panic). All these ef- 
fects are simulated by varying of the model parameters. 
Model obtained saves this opportunities. Modifications and 
improvements made here mainly concern the quality of 
pedestrians behavior reproducing. They allow more realis- 
tically (carefully) simulate analysis that people accomplish 
while they choose a direction for moving. Idea of parameters 
adaptivity makes model more flexible and closer to real life. 

At first let us consider components that determine proba- 
bility pij. In contrast to FF model one can distinguish two 
different types of terms in (Q]i. 

Term Ajj characterizes the physical possibility to move, 
Ajj E [0,1]. It takes maximal value if movement conditions 
in the direction are favorable. And Ajj = if there is no free 
space to move. Term Ajj proportionally decreases with the ad- 
vent and approaching of some obstacles (people, wall, etc.) in 
the direction. 

Other terms e ksSi ' , e koD 'j, e k ' vary form 1 to ™ (in gen- 
eral case) and characterize style of people behavior. Mini- 
mal value of parameter (k$ = 0, ko = or k\ = 0) means 
that correspondent feature of behavior isn't realized and term 
doesn't affect the probability. If all three terms are minimal 
then pedestrians walk free. And in this case only term Ay de- 
termines the transition probability for each next-neighbor cell 
in accordance with people features: keeping apart from 
other people and obstacles, patience. 

In FF model fl][3][5 1 pedestrians stay at present cell if there 
is no space to move only. Here we give pedestrians the oppor- 
tunity to wait when preferable direction will free even if other 
directions are available for moving at this time. Such behav- 
ior is reside to low and middle densities. To realize it tran- 
sition probabilities ([]} don't include a checking if cell 
occupied or not. In this case transition probabilities (fTJ can 
be considered as a rate of wish to go to the certain directions. 
Possibility not to leave present cell is realized in step 3. 

The idea of spatial adapted parameters is introduced here. 
It's clear that conditions can' t be equal for all people involved. 
Position of pedestrian in the space may determine some fea- 
tures of the behavior. Thus clogging situations are natural and 
well pronounced for nearest to exit (destination point) space. 

For other areas it's not typical but it's possible. To realize it 

s 

the coefficient m J i ' s . for clogging parameter fj is introduced. 

SIMULATION RESULTS 

In order to test our model a regular evacuation process was 
simulated. This means that ko — in examples presented. 
And we set kj = 0, /j = 0. 

One pedestrian case 

There was simulated evacuation of one person (n=l) from 
a room 6.8m x 6.8m (17 cells x 17 cells) with one exit (0.8m) 



in the middle of a wall. Recall that the space is sampled into 
cells of size 40cm x 40cm which can either be empty or oc- 
cupied by one pedestrian only. Static field S was calculated 
in accordance with (3J. Stating position is a cell in a corner 
near wall opposite to the exit. Pedestrian moves towards the 
exit with v = v mux = 1 . For such sampled space minimal value 
of time steps that require to leave the room starting from ini- 
tial position is r m ,„ = 26. 

Different combinations of parameters k$ and r were consid- 
ered. Total evacuation time and trajectories were investigated. 
Following table contains results over 500 experiments. 

Table 1. Modes T mo of total evacuation time distributions in 
one pedestrian case. 
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Figures[T]show total evacuation time distributions for some 
couples of the parameters from table[T]over 500 realizations. 
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e)k s = 4, r=l. f) k s = 4, r = 17. 

Figure 1. Total evacuation time distribution for different k$ 
and r over 500 experiments. 
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e)k s = 4,r=l. f)*s = 4,r=17. 

Figure 2. Tracks of one pedestrian for different and r over 
500 experiments. 



the direction to destination point. In other words one can say 
that pedestrian doesn't see, knows, and wants not so much to 
go to destination point (exit). Last one (r = 17, ks = 4) de- 
scribes situation when pedestrian sees, knows, and wants to 
go to destination point very much. 

Note that if r — 1 model presented corresponds to FF- 
model [ 1 3 1 with the same other parameters. And it's clear 
that the pedestrian patience that was introduced in the model 
doesn't pronounced in one pedestrian case. 

One can see that for small ks = 1 mode T mo is very depen- 
dent on parameter r. The bigger parameter ks is an influence 
of r to T mo is less pronounced. But the bigger parameter r 
is more natural way pedestrian chooses — tracks are more 
close to a line connecting starting point and exit (non the less 
random component takes place). Under r > 1 a proximity of 
wall decreases the probability ([TJ to move in this direction. 
Thereby tracks are forced to tend to natural one. Thus param- 
eter r fulfils its role to simulate environment analysis here. 



Figures |2]present tracks of pedestrian ways over 500 real- 
izations for some couples of the parameters from table Q] 

Different moving conditions are reproduced by these com- 
binations. They vary from r = 1, ks = 1 to r = 17, ks = 4. 
The former case can be interpreted as pedestrian moves in a 
very low visibility (or by touch) but he approximately knows 



Many pedestrians cases 
Influence of kg and r 

For this collective experiment the space was a room 16m x 
16m (40 cells x 40 cells) with one exit (0.8m). Initial num- 
ber of people is = 300 (density p « 0.19). Initial positions 
are random and people start to move towards the exit with 
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a) k s = l,r= l,T total =509 
t=20 f = 210 
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The problem comes from computational aspect of coefficient 
Aij. It's a positive that in the model obtained people avoid to 
approach to walls. But this effect has to be less pronounced 
with approaching to wall (walls as in a subsection example 
below) surrounding exit. Thus parameter r has to be spatial 
adaptive. 

Next table demonstrates numerical description of cases 
presented. Let f a be frequency to choose direction a = 
{N,E ,S,W,C} over all experiment for each couple of param- 
eters. Here N,E ,S,W,C are north, east, sought, west, center 
(stay at present place) correspondingly, M — total number 
of movements including stayings at current position over all 
experiment. 

Table 2. Direction frequencies. 
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a) k s = 3, r = 1, T total = 336 
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bH 5 = 3,r = 40,r,„, o/ =317 
Figure 4. Evacuation stages for 300 people for ks = 3 and 
different r. 



v = Vmux = 1 ■ Exit is in the middle of east wall. Figures [3] |4] 
present typical stages of evacuation process for different ks 
and r. 

One can see that evacuation dynamics in case a) differs 
from case b) in both figures [3] [4] The reason of it is different 
parameters r. If r > 1 people avoid to approach to walls. And 
while crowd density allows pedestrians try to follow more 
natural way to the exit. The bigger r is shape of crowd in 
front of exit is more diverse from the case of r = I. (Note if 
r = 1 we have FF-model with the same other parameters.) 

One can notice that the closer to the exit pedestrians are 
circle shaped crowd in front of the exit is more unrealistic. 



Comparing cases 1) and 2), 3) and 4) correspondingly one 
can notice that greater r leads to significant redistribution of 
flow. In the case of r = 1 low fc says that people move as 
much as possible (because people can stay at current position 
if all nearest cells are occupied only). If r = 40 the opportu- 
nity to wait is realized — in cases 2) and 4) fc has the greatest 
value, fw has the smallest value (this direction is opposite to 
the exit), fy, fs are approximately equal because exit is in the 
middle of the wall. Thus in cases 2) and 4) model reproduces 
more natural decision-making process. 

At the same time let us remark that increasing ks makes 
evacuation process more directed and reduces total evacua- 
tion time (see cases 1) and 3), 2) and 4) correspondingly). 

Influence of exit position 

For this collective experiment the space was a room 6.8m x 
11.2m (17 cells x 28 cells) with one exit (0.8m). Initial 
A' = 150 (density p w 0.31). Initial positions are random and 
people start to move towards the exit with v = v max = 1. 
Spaces are presented in the figures [5h,[5]3. 

Two combinations of parameters ks and r were considered. 
And total evacuation time was investigated. Table [3] contains 
results over 100 experiments. Initial positions of people were 
the same for all experiments. 

One can notice once again that in the case a) parameter r 
doesn't influence on total evacuation time under such ks- But 
in other case increase of r leads to significant delay of evacu- 
ation. The reason of it is "computational" repulsion from the 
walls. And as a result pedestrians don't use corner between 
wall and exit, exit is not fully used, total evacuation time in- 
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a) b) 
Figure 5. Initial people position and exits position. 

Table 3. Modes T mo of total evacuation time distributions 
with middle exit (a) and corner exit (b). 
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creases. So this is one more example that shows necessity in 
at least parameter r adaptivity. 

CONCLUSION AND FURTHER PLANS 

In the paper the intelligent FF cellular automation model 
is presented. Modifications made are to improve realism 
of the individual pedestrian movement simulation. The fol- 
lowing features of people behavior are introduced: keeping 
apart from other people (and obstacle), patience. Idea of spa- 
tial adaptation of model parameters is pronounced and one 
method is presented (parameter £ty). Model obtained saved 
opportunities to reproduce variety of collective effects of 
pedestrian movement from free walk to escape panic and took 
more flexibility. 

The simulation made showed real improvements in 
decision-making process in comparison with basic FF model 
and pointed out some problems. The following points seem to 
be very important for realistic pedestrian simulation and are 
under future investigations. 

Jam 

In a case of emergency appearing of clogging situation in 
front of exit often leads to appearing fallen or injured people 
(or jam). The physical interactions in such crowd add up and 
cause dangerous pressures up to 4,450./Vm~ 1 [2] which "can 
bend steel barriers or push down brick walls". Fallen or in- 
jured people act as "obstacles", and escape is further slowed. 
Continuous model can reproduce pushing and physical inter- 
actions among pedestrians. CA model doesn't allow to do it. 
Parameter ftjj works in the model as some kind of local pres- 



sure between the pedestrians (the higher fijj is pedestrians are 
more handicapped by others trying to reach the same target 
cell). But by means of fallen or injured people are not 
simulated. Our further intention is to produce method to eval- 
uate common pressure to each pedestrian in CA model, and 
if pressure is over some barrier to indicate correspondent cell 
as new obstacle. 

Parameters adaptation 

There are at least two reasons for parameters to be adaptive. 
One of them is learning that is reside to people. Therefor at 
least parameters ks and kg need to be time adaptive and spa- 
tial dependent as well. Parameter r needs to be spatial adap- 
tive because of negative computational effects. So methods to 
adapt model parameters are under further investigation. 
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